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A relativistic neutral scalar field is investigated on the basis of the Schwinger-Dyson 
equation in the non-equilibrium thermo field dynamics. A time evolution equation for a 
distribution function is obtained from a diagonalization condition for the Schwinger-Dyson 
equation. An explicit expression of the time evolution equation is calculated in the \<fi 4 
interaction model at the 2-loop level. The Boltzmann equation is derived for the relativistic 
scalar field. We set a simple initial condition and numerically solve the Boltzmann equation 
and show the time evolution of the distribution function and the relaxation time. 

§1. Introduction 

The relativity and the quantum field theory provide a general basis for under- 
standing a model of elementary particles at high energy. The theory should be ex- 
tended to include a statistical thermodynamics for macroscopic phenomena of many 
particles systems. Some formalism for thermal quantum field theories is proposed 
to evaluate in and out of equilibrium systems. The real-time formalism is necessary 
to include time evolution of the system. Y. Takahashi and H. Umezawa proposed 
the thermo field dynamics (TFD) which is the real-time formalism based on the 
canonical quantization.^ - ® After that TFD has been applied to various physical 
systems. 

In 1985 T. Arimitsu and H. Umezawa extended TFD to a non-equilibrium sys- 
tem with spatially homogeneous distributions. 4 '" 5 ^ It has been found that the time 
development of a particle number density is derived by the self-consistent renormal- 
ization condition^'®"® An equation with the structure of the Boltzmann equation 
is derived in the non-equilibrium TFD (NETFD). An alternative procedure to derive 
a Boltzmann-like equation has been proposed by evaluating the expectation value of 
the number operator in terms of a renormalized field in RefsP'^ 1 It is also obtained 
by diagonalizing the Green's function at the equal time limit instead of the thermal 
self-energy diagonalization scheme in Ref! 1 - 1 ^ TFD for spatially inhomogeneous non- 
equilibrium system has been studied in RefsP^ NETFD successfully describes 
macroscopic phenomena in a low energy non-relativistic system. 

It is also quite interesting to make investigations various phenomena associated 
with non-equilibrium dynamics at high energy. As far as we know, only a little 
work has been reported in the study of NETFD for a high energy relativistic system. 
P. Henning introduced independent weighting functions for relativistic fields with 
positive and negative energy and described an extended Bogoliubov transformation 
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for thermal doublets; m NETFD was applied for a relativistic scalar field with a 
four-point self-interaction in RefPSJ A renormalization condition is imposed on the 
thermal self-energy. Thus a Boltzmann-like equation is derived from the expression 
for the particle number density of the scalar field. 

In the present paper we would like to find a natural derivation of the self- 
consistency conditions available for relativistic scalar fields by using the Schwinger- 
Dyson (SD) equation. In Sec. 2 we briefly review NETFD for a scalar field and 
introduce the thermal Bogoliubov matrices whose parameter defines a particle num- 
ber density. In Sec. 3 we evaluate the SD equation in NETFD. A Boltzmann equation 
is derived from a time evolution of a parameter in the thermal Bogoliubov matri- 
ces. In Sec. 4 a scalar field with a Xcp interaction is investigated according to the 
NETFD. We calculate the thermal self-energy at 2-loop level. In Sec. 5 we give an 
explicit expression for the Boltzmann equation in the Xcf> interaction model. Solving 
the Boltzmann equation starting from a simple initial state numerically, we show the 
time evolution of the particle number density and the relaxation time. We discuss 
how the off-shell mode contributes the relaxation process. Some concluding remarks 
are given in Sec. 6. 

§2. Non-equilibrium thermo field dynamics for relativistic fields 

There are several formalism to introduce the thermal dynamics into the quantum 
field theory. In TFD the statistical average is replaced by an expectation value in a 
pure state called the thermal vacuum. The thermal vacuum is defined by extending 
the Fock space structure. A creation and an annihilation operators, a\ a are doubled 
by introducing tilde operators, a\ a. Thus the commutation relations for the bosonic 
creation and annihilation operators are extended to be 

[a p ,a q ] = (2nf5^(p-q), (2-1) 

[dp,^ q ] = (27rf5^(p-q), (2-2) 
others = 0. (2-3) 

The time evolution of the non-tilde operators are generated by an ordinary non- 
thermal Hamiltonian, H, for the considering system. In a similar manner the time 
evolution for the tilde operators is described by the tilde conjugate Hamiltonian, H. 
In defining the tilde-Hamiltonian, H, we use the following tilde conjugation rules, 

(AtA 2 r = A X A 2 , (2-4) 
( Cl Ai + c 2 A 2 r = c\Ai + c* 2 A 2 , (2-5) 
(Ay = A, (2-6) 

where c\ and c 2 are c-numbers, and A± and A 2 are arbitrary operators. The tilde- 
Hamiltonian, H, is constructed by only the tilde operators, a and a) . The time 
evolution of both the non-tilde and tilde operators is described by the total Hamil- 
tonian of the system, 

H = H — H. (2-7) 
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The Fock space is extended to the state space spanned by both the non-tilde 
and tilde creation operators. The thermal vacuum is defined by the thermal Bogoli- 
ubov transformation of the eigenstates of the Hamiltonian. Non-equilibrium degree 
of freedom can be also introduced in TFD through the thermal Bogoliubov trans- 
formationPS® 

^-^■ ( = BK(t)) tt "<w, 

e«e^-' = a£(i)iT VpW)^, (2-8) 



where uj p = \/p 2 + m 2 is the relativistic energy eigenvalue for bosonic fields with a 
momentum, p, and mass, m. The upper indices are defined by the thermal doublets 
notation, 

a p = ( ~t J ' S p = ( a P ~"-p I ' ( 2 ' 9 ) 



It is assumed that the thermal Bogoliubov matrices, B{n p {t)) a ^ and B~ l (n p (t)) a " 
have the same forms with the ones in equilibrium, 

%(!))■» = ( 1 ~ n f ] 

r 'W')r' = (l (2 ' n) 

The Bogoliubov parameter, n p (t), depends on time and the relativistic energy eigen- 
value through the momentum, p. 

In NETFD the canonical quantization has not been fully established for rela- 
tivistic fields yet although some attempts have been doneP^J In the present paper it 
is assumed that the quantum correction is calculated by perturbative theory in the 
Fock space spanned by the ^-oscillators. The thermal vacuum, \9), is defined by the 
transformed operator, £ p , 

<0|4 = 6|0>=O, (0\it = i P \0)=0. (2-12) 
The state, \9), is invariant under the tilde conjugation (|2-4j) . (|2-5p and (12 -6p . 

i(e\r = (oi mr = \o). (2-13) 

In equilibrium system the statistical thermal average is obtained by evaluating the 
expectation value under the thermal vacuum. 

A neutral scalar field is expanded by using the a-oscillators, i.e. inverse of the 
thermal Bogoliubov transformations (|2-8p on the ^-oscillators, 

m = I pi^^) eip ' s + a l^~ ip ' x i ( 2 - 14 ) 

~ m = I + 4(t*y p ' x )- ( 2 - 15 ) 
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From Eq. (|2-7p the time dependence for the operators a and a) is identical to a) and 
a, respectively. Thus the thermal doublets notation is defined for the neutral scalar 
fields by 

d 3 P 1 / f ap(t x ) \ Jp . x f a P {t x ) \ 



= / " " ~ U J a e^ x + e 
' (2vr) 3 y2Z^ 1 V 4(tx) J V a p (t x ) 1 



d 3 p 



(2vr) 3 y2^ 
d 3 p 



{a%(t x yP x + {r,a p {t x ) T r^P x } , (2-16) 



r(x) = J ( ^ 3 ^1- { (4(i,) - gpfe)) + (a p {t x ) - al(t x )) jP x ) 

d 3 p 



(2^)3^/2^ 



{a^Qe-V-* + (a p (^)r 3 ) Q! e i ? , - a; } , (2-17) 



v 



where T3 is the third Pauli matrix acting on the thermal doublets. 

Below we express fields and propagators by the t-representation, which is defined 
by taking the spatial Fourier transform. A thermal propagator is defined by an 
expectation value of the time-ordered product of two fields in the thermal vacuum. 
It has a 2 x 2 matrix form with respect to the thermal index. For a neutral scalar 
field it is defined by 

D a \t x - t y ;p) = {9\T[r(t x ;p)^(t y ;p)]\6), (2-18) 

where a and ft are the thermal indices. Thus the thermal propagator is obtained by 
the thermal Bogoliubov transformation of the non-thermal onep^^'^'® 

B a \t x - t y ;p) = B-\n p {t x )) a ^DlJ(t x - t y ;p)B(n p (t y ))^ 

+{T 3 B(n p (t x )) T }^D^(t x - t y -p){B-\n p {t y )) T T^, (2-19) 

where -Dj 1 ^ 2 ^. — t y ;p) and D^]^{t x — t y ;p) are the retarded and the advanced 
parts for the non-thermal propagator which are given by 2 x 2 matrices with thermal 
indices 71 and 72. The non-thermal propagators for the scalar field have diagonal 
forms. Each component in Dq ^ and Dq a is given by 



Doln(t x 


-t y ;p) 


= 0(t x - 




-ty]p) 


= -Gity 


Dl] A (t x 




= 0(ty - 


Do 2 A (t x 


~ty-,p) 


= -e{t x 



j_ -iU)p-{t X -ty) 



2u p 



(2-20) 



1 Q-iUpitx-ty) ^ (2-21) 



2ui p 
2u> p 



(2-22) 



1 e iuj p -(t x -t y ) ^ (2-23) 



2u p 



other components = 0. 



In NETFD we can choose the boundary condition for the thermal Bogoliubov 
matrices (|2-11|) to derive the ordinary Feynman rules for relativistic scalar fields 
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It is convenient for perturbative calculations. In TFD the thermal average of the 
particle number density is given by the expectation value of the number operator 
a} R p an,p in the thermal vacuum, 

n p (t) = (6 H \a^ p (t)a HtP (t)\6 H ), (2-24) 

where the fields and ground state with the lower index, H, are written in the Heisen- 
berg representation. It coincides with the Bogoliubov parameter, n p , in the Bogoli- 
ubov matrices (|2-11|) in an equilibrium state at tree level. Since quantum corrections 
can induce instability of the thermal vacuum (|2T2p . the Bogoliubov parameter does 
not always correspond to the particle number density out of equilibrium states. Be- 
low we assume that a stable vacuum is found by redefining the Bogoliubov parameter 
and the observed particle number density is obtained by the Bogoliubov parameter 
in the stable vacuum at the equal time limit P 



§3. Time evolution equation 



H. Umezawa and Y. Yamanaka have introduced the self-consistent renormal- 
ization condition to derive the Boltzmann equation in NETFD P The hat- 
Hamiltonian (|2-7p does not contain terms proportional to ££. However, such a term 
is induced through the quantum correction in the thermal self-energy out of equi- 
librium states. It breaks the condition (|2-12p for the thermal vacuum. The terms 
proportional to can be eliminated by the thermal Bogoliubov transformation. 
We impose that the terms proportional to vanish at the equal time limit. The 
perturbed thermal Bogoliubov parameter can be fixed by this self-consistent renor- 
malization condition. The Boltzmann equation appears as a consequence of this 
renormalization condition. It is also derived from diagonalizing the full propagator 
given by the Dyson equation at the equal time limit in non-relativistic quantum field 
theories.^ 1 In this work we improve this diagonalization condition to be suitable for 
a relativistic scalar field. 

In TFD the SD equation for the scalar field is given by 

D°^{t x -t y yp) = D^(t x -t y ;p) 
+ / dt Zl dt Z2 D^(t x - t Zl ;p)i£^(t Zl - t Z2 ;p)D^(t Z2 - t y ;p), (3-1) 

where D^f (t x —t y ;p) denotes the full thermal propagator, D^{t x — t y ;p) the thermal 
propagator at tree level and (t zi — t Z2 ; p) the full thermal self-energy. Another 
expression for this equation is 

Dff(t x -t y ;p) = D^(t x -t y ;p) 
+ J dt zl dt Z2 D^(t x - t zl ;p)iS^(t Zl - t Z2 ;p)D^(t Z2 - t y ;p). (3-2) 

We suppose that the full propagator can be decomposed in a similar form with the 
tree level one (|2-19p . 

D^(t x - t y ;p) = B-\n HtP {t x )) a ^D^(t x - t y ;p)B(n H , P (t y ))^ 
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+{T3B(n H , p (t x )) T }^Dl}J(t x - t y] p){B-\n H , p {t y )) T T^i 



(3-3) 



where nn,p{t x ) an d n>H, P (ty) are the Bogoliubov parameters acting on the full thermal 
propagator from the left- and the right-hand sides, respectively. For simplicity we 
omit the thermal indices and the momentum label p in the thermal propagator and 
the Bogoliubov parameters below. 

We decompose the propagators in the SD equation (|3-ip in accordance with 
Eqs. (|2-19p and (|3 3[) to see the matrix structure of the thermal propagator. After 
the Klein-Gordon operator, {d 2 x + p 2 + m 2 ), is applied on the left to both sides, 
Eq. §FQ) reads 

(8l +p 2 + m 2 ) [B-^nnit^DH^ - t y )B{n H {t y )) 

+T 3 B{n H (t x )) T D HA (t x - t y )B- 1 {n H {t y )) T T ?J 
= (dt +P 2 + m 2 ) [B~\n(t x ))Do, R (t x - t y )B{n{t y )) 

+T 3 B(n(t x )) T D 0A {t x - t y )B- 1 {n{t y )) T T 3 

+ J dt Zl dt Z2 {d 2 tx +p 2 + m 2 )\B- 1 {n(t x ))D ,B.(t x - t Zl )B(n(t Zl )) 

+T 3 B{n{t x )) T D 0A {t x -t Zl )B- 1 {n{t Zl )) T T^iE{t Zl -t Z2 ) 
x B- l (nH{t Z2 ))DH,R{t Z2 -t v )B{nH{ty)) 



+T 3 B(n H (t Z2 )y D HA (t Z2 - t y )B- L (n H (ty)Y r 3 



(3-4) 



where n(t) is the Bogoliubov parameter for the non-perturbed operator included in 
the thermal propagator at tree level. 

The Bogoliubov matrices and the Klein-Gordon operator do not commute. Sub- 
stituting the expression (|2-11|) to the Bogoliubov matrices and inserting an identity 
operators, B~ 1 (nH)B(riH), we obtain 



B {n H {t x )){d tx +p +m)[ Q 



,propl (t x i ty , 

+r 3 B(n H (t x )ndl+p 2 + m 2 )( °»> A %~ V D 22 , ° , 

\ ^ A,propl\ L xi b y ) ^ H A\ x V) 



- / dt. 



B-\n H {t x )) 



2^R.(tx t s )D^ft(t s ty) 9xl{t x .,t y ,t s ) 
2^A(tx ~ t s )D 2 ^ ft(t s ~ t y ) 



+B-\n H (t x )) 

+T 3 B(n H (t x )) T 

+T 3 B(n H (t x )) T 



9x2 {t x ,ty,tgj 



\^h{tx ~ ts)Dff j{(t s — ty) 
2^A{t X —t S )D Hj ^(t s —ty) 

— ^YtA^tx ~ ts)D 2 ^ft(t s — ty) 

~2~2jR{tx ~ ts)D^jft(t s — ty) gx3(tx,ty,t s ) 





B(n H (t y )) 
B(n H (t y )) T T 3 
B(n H (t y )) 
B-\n H {t y )) T T 3 



B(n H (t y )) 



2^A(tx ts)D]}^(t s ty) „ . )) T T 3 

9x4{tx,ty,t s ) l^Sjl{t X — t S )D j^j^itg — ty) J ^ 
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-l6(t x - ty), 



(3-5) 



where the thermal self-energies for the retarded and the advanced propagators, E R 
and S A , are defined by 



U R = U 11 + U 12 = U 21 + U 22 , 



E A = U u - U 21 = S 22 - u 12 . 



(3-6) 



The off-diagonal elements, OrmWojjIj m the propagator (|3-5p are defined to satisfy 
the following equations, 



(d 2 x + p 2 + m 2 )0 R , prop i{t x ,ty) 

= h H (t X )D% R (t X - ty) + 2h H {t X )(d tm D% >R {t X - ty)) 
-h(t X )D 22 R (t X - ty) - 2h(t X )(d tx D 22 R (t X - ty)), 

(d 2 x + p 2 + m 2 )0 A , P ropi(t x ,ty) 

= -n H (t x )D% A {t x - t y ) - 2h H (t x )(d tx D 1 H 1 A (t x - ty)) 
+n{t x )Dl] A (t x - t y ) + 2h(t x ){d tx Dl] A {t x - ty)), 

We define the off-diagonal elements, g x \ ~ g x ±, by 

9x1 (fx lty,t s ) 

= ^^{tx ~t S ) + h-(t X ,t S )^D 2 £ iR (t S - ty) 



dt. 



{h(t x )D 22 R (t x - t s ) + 2h(t x )(d tx D 22 R (t x - t.))} 



x iU A (t s - t z )Djf iR (t 8 - t 

9x2 \t x ityi^s) 

1 r 



2 {^ n (** - t a ) + h+(t x ,t s )}D% ;A (t s - ty) 

{h(t x )D 22 R (t x - t s ) + 2h(t x )(d tx D 22 R (t x - t t ))} 



x iS A {t s -t z )D^ A {t z -t y 
1 r 



^ 22 (.tx ~t S )+ h + (t X ,t S )}D 2 ^ R (t S - ty 



dt. 



{n{t x )Dl] A {t x - t 3 ) + 2h{t x ){d tx Dl] A {t x - t a ))} 
x iU A (t s — t z )D 2 ^ R (t s — ty) 



(3-7) 



(3-8) 



(3-9) 



(3-10) 



(3-11) 
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9xi(t x , ty,t s ) 

= ~^{Z 21 (t X - t S ) + h^Q^D^ts - ty) 



dt } 



n.( 



with 



{t x )Dl] A {t x - t s ) + 2ii{t x ){d tx Dl] A {t x - 
x iU A (t s ~ t z )D l H A {t s - ty) 



h-(t, t') = n H (t')S R (t - t') - n H (t)S A (t - if), 
h + (t,t') = n H (t')U R (t - t') + n H (t)E A {t - t'). 



(3-12) 



(3-13) 
(3-14) 



It should be noticed that the first and the second derivatives of the Bogoliubov 
parameters appear in Eq. (|3-5|) through the off-diagonal elements, R r A \ prop i, and 
9x1 ~ 9x4- 

We also rewrite the SD equation (|3-2p . Applying the Klein-Gordon operator 
from right to both sides of Eq. (|3-2p . we obtain 



+T 3 B(n H (t x )) T D HA (t x - t y )B- 1 (n H (ty)) T T 3 {tl + p 2 + m 2 ' 



B-Hnnit^DH^ - t y )B(n R (t y )) 



B- L (n(t x ))D 0tR {t x - t y )B(n(t y )) 



+r 3 B(n(t x )) T D 0A (t x - ty)B- L (n(t y )) T T 3 (d 2 +p I + m I ) 



+ I dt z , dt 



B-\n H (t x ))D H , R (t x - t Zl )B(n H (t Zl )) 



+T 3 B(n H {t x )) T D HA {t x - ^ 1 )S- 1 (n / f(t 2l )) T r 3 
xiU(t Zl — t Z2 )\B-\n{t Z2 ))D , R {t Z2 - t y )B(n(ty)) 
+T 3 B(n(t Z2 )) T D 0A (t Z2 -t y )B- 1 (n(t y )) T T 3 \( i d 2 ty + p 2 + m 2 
Substituting Eq. (|2-lip to this equation, Eq. ()3-15p reads 

Bjj R {t x ty) Ri p r0 p2(tx,ty) 



(3-15) 



B-\n H (t x )) 
+T 3 B(n H (t x )) T 

dt 



{di+p 2 + m 2 )B(n H {t y )) 



D 22 R {t x -t y ) 

D HjA (t x — t y ) ^ \ _j_ m 2\_g-l^ nfl .^ )) T T3 

OA,prop2{txjty) Bfj A (t x — ty) J * s 

2^^H,R^ X ~ ^s)^r{Ps ~ ty) 9yl(t x , ty, t s ) 
^Bfj R {tx ~ t s )^A{ts ~ ty) 



B-\n H (t x ))[ ^""' () ' "T^i'u^^ZTi )B(n H (t y )) 



9y2(tx, ty, t s ) ~ ^Djj R {t x — t s )U R (t s — ty) 



^ ' ' ' R^ x t s )U A (t s ty 







B-\n H (t y )) T T 3 
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+ 3 ( h( x)) — ^D 2 ^ A (t x — t s )S R {t s — t y ) g y 3{t x ,t y ,t s ) ) ^ h( y)) 



+r 3 B(n H (t x )) q 



~ ts)^A{ts ~ ty) 
9yi{t x , ty, t s ) 2^H,A^fx ~ t s )Z2 R (t s — ty 



B-\n H {t y )) T T 3 



-iS(t 



X ty) 



where the off-diagonal elements, O r ^ , P rop2> are denned to satisfy 
Rt p r0 p2{tx,ty)(d 2 ty + p 2 + m 2 ) 

= -h H {ty)D l l R (t X - ty) - 2h H (ty)(d ty D% R (t X - ty)) 

+n(t y )Dl] R {t x - t y ) + 2n{t y )(d ty Dl] R {t x - t y )), 
A ,prop2(t x ,t y )( d 2 ty + p 2 + m 2 ) 

= n H {ty)D 22 A (t x - ty) + 2h H (ty)(d ty D 22 A (t X - ty)) 

-Htx)D 22 A (t x - t y ) - 2h(t x )(d ty D 22 A (t x - ty)), 
and the off-diagonal elements, g y \ ~ g y ±, are 
9yi{t x , t y ,t s ) 

= \D l i R (t X - t s )[S 12 (t s - ty) + h-(t 8 ,ty)} 

dtz B>H R (t x - t z )iU R (t z - t s ) 



x{n(ty)D£ R (t S - ty) + 2n{ty){d ty Dl] R {t S - ty))} 

Qy2^t X , t y , t s ) 

= ^iRi^ - ts){Z U (t s - ty) + h+(t S ,ty)} 



f 

J — l 



dt. 



DH, R (t x - t z )iU R (t z - t s ) 



x{h(ty)D 22 A (t s ~ ty) + 2h{ty)(d ty D 22 A (t S ~ ty))} 

9y3(t x , t y ,t s ) 

= ~\D 22 A (t X - ts){S 22 (t S - ty) + h + (t S ,ty)} 

D 2 £At x -t z )iE R {t z -t s ) 



"OO 

+ I dt z 

-oo 



{Hty)D^ R (t s - ty) + 2h{t y ){d ty Dl] R {t s - t y ))} 



(3-16) 



(3-17) 



(3-18) 



(3-19) 



(3-20) 



(3-21) 
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9yi(txi ty,t s ) 

-D 2 £ A (t x - t s ){£ 21 (t s - ty) + h„(t s ,ty)} 



2 

+ / dt 



\n{ty)Dl 2 A {t S - ty) + 2h(ty)(d ty D 22 A (t S - ty))} 



(3-22) 



For non-relativistic fields the particle number conservation law is derived from 
the difference between the SD equation applying the differential operator, idt x + 
p 2 /2m, on the left and right sides at the equal time limit We adapt the procedure 
to relativistic scalar fields. Subtracting the Eq. (|3-16p from Eq. (|3 5|) . we obtain 



B-\n H {t x )) 



( af +p 2 + m2)/ D% tR (t x t yJ 







Ofl^propl (t x , ty) 



^H,r(.^x t y ) 



Ojl j p r op2(txi ty) 
Djinitx — ty) 



(d 2 t +p 2 + m 2 



dt. 



2^R(tx tg)DjJ R (tg ty) 9xl(t x , ty, t s j 







h^A{t x — t s )Dj^ n(t s — t 



+r 3 B(n H (t x )) 



R(t X t S )Sjl(t S ty) Qyl(t X -,ty,t S/ 

R^P X ~ t s )^A(ts ~ ty) 



T 



(Of +p 2 + m 



B(n H {t y )) 




O A,propl(t x ,ty) Dff J^{t x t 



O A,prop2{t x ,ty) D 2 ^ A (t x ty) J ^ *» 



dt 



2^A(t x t s )D^j A (t s t 







9xi(txity,t S ) 2~UR(t X tg)Djj A (tg ty) 

A^P X ~ ts)^A(ts ~ ty) 
9yi{t x , ty, t s ) ^HA^ 1 — ts)^R(ts ~ ty t 

9x2(t X ,ty, t S ) — ^Uji(t X — t s )D 2 ft A (t S — ty) 

-\E A {t x -t s )D^ A {t s -t y ) ' 



B-\n H {t y )) T Tz 



-B- l {n H {t x )) / dt 



9y2 (t x ,ty,t s ) 2 
~\^H,R^ X ~ ts)^*A(ts ~ ty 



1 R^x — t s )^R(t s — t y ) 



B- L (n H (ty)YT 3 



-r 3 B(n H (t x )y / dt 



— ^UA(t X — t S )D 2 ^ R (t S — ty) 

'2^R{tx ~ tsjDjj p>(t s — ty) 9x3{txi ty, t s ) 
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I^ffife — ts)^A(ts ~ ty) 



1 -^li.A^x ~ t s )2J R (t s ~ ty) 9y3(txity,t s 



B{n H (t y )) = 0. (3-23) 



Each elements of the matrices between the Bogoliubov matrices t^B t {- ■ ■ )B and 
B~ l (- ■ ■ )B~ 1T Tz and the diagonal element in B~ 1 (- ■ ■ )B and TsB T (- ■ ■ )B~ 1t t^ sat- 
isfies trivial equation, = 0, at the equal time limit, t x — > t y , in the case of the 
relativistic scalar field with a four-point self-interaction, as is shown in the next 
section. The remaining off-diagonal elements should satisfy 



lim 



(d 2 x + p 2 + m 2 )0 R:PTOpl (t x ,t y ) - RjProp2 (t x ,t y )( d \ y + p 2 + 



m 



dt, 



^9x1 (j^X J ty, tg) (jy\ (t X , ty , t S ) J' 



lim 



(df x + p 2 + m 2 )OA,propl(tx,ty) - A ,prop2(tx,ty)( 8 \ y + p 2 + 



m 



- / dt. 



^ 9xi(tx) ty, t s ) g y i(tx,t y) t s )} 



0. 



(3-24) 



(3-25) 



Substituting Eqs. (ETSl) . (f3T71> and (|3~T9]l into Eq. (^241 . we obtain 



lim [h H (t x )D HR (t x - t y ) - h(t x )D 0R (t x - t y 



+ n H (ty)D H R (t X - ty) ~ rl(ty)D R (t X - ty) 

+2n H (t x )(d tx D 2 l R (t x - t y )) - 2h(t x )(8 tx D 22 R (t x - ty)) 

+ 2h H (ty)(d ty D 1 H \ R (t X - ty)) ~ 2n(ty)(d ty Dl] R (t X - ty)) 
^2 j dt s[{^ l2 (tx-tg)+ih^(t X ,tg)^D 2 ^ R (t S -ty) 

-Dn )R (t x - t s ){iZ 12 (tg - t y ) + ih„(t s ,ty)} 



0. 



(3-26) 



This equation describes the time evolution of the Bogoliubov parameter, n(t). 

We write the differences between the perturbed and the unperturbed Bogoliubov 
parameter as 



v(t) = n H (t) - n(t). 



(3-27) 



Substituting Eqs. (fA^ - fjAdHI) to Eq. (^26\\ . a delta function , 5(t x - t y ), ap- 
pears from the first derivative of the propagator. It diverges at the equal time limit. 
The divergence is canceled out in Eq. (|3-26|) . The ordinary Boltzmann equation does 
not contain terms proportional to the second derivative of the thermal Bogoliubov 
parameter, v. These terms in Eq. (|3-26[) are canceled out at the equal time limit. 
Therefore the equation (|3 26|) is simplified to a Markovian equation, 



>(t x ) = ~Yim t j dt s [{i£ 12 (t x - t s ) + ih-(t x ,t s )}D% iR (t 8 - t y ) 

-DH :R (t X - tg){iZ 12 (tg - ty) + ih-(t a ,ty)} 



(3-28) 
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In a similar manner Eq. ()3-25[) reduces to a Markovian equation for the advanced 
propagator, 

H*x) = ~ t jim J dt s [{i£ 21 (t x -t s )+ih_(t x ,t s )}D l ^ A (t s -t y ) 

-D 2 £ A (t x - t s ){i£ 21 {t s - t y ) + ih-(t„ty)}] . (3-29) 

Therefore the time evolution of the Bogoliubov parameter is determined by solving 
Eqs. (f3^28l) and (^291) . 

§4. Boltzmann equation in a A0 4 theory 

It is expected that the time evolution of the number distribution is given by 
the Boltzmann equation in NETFD. In this section we show that the Boltzmann 
equation is obtained from Eqs. (|3-28|) and ()3-29p . We consider a neutral scalar field 
with a four-point self-interaction in (1+2) dimensions, for simplicity. To derive the 
Boltzmann equation we calculate the self-energy and evaluate the time evolution 
equations (|3-28p and ()3-29p . We start from the Lagrangian density 

= \(d,cp{x))(d^{x)) ~ \m 2 m 2 ~ j { 4>{xf. 

(4-1) 

In TFD the total Lagrangian density, C, for non-tilde and tilde fields is given by 

t{x) = C(x) - C(x), (4-2) 

where C is the tilde conjugate of the Lagrangian density C. In this model the self- 
energy at 1-loop level can not contribute to the right-hand sides in Eqs. (|3-28|) and 
(|3-29p . We have to calculate the self-energy at the 2- loop level. 

The self-energy is calculated by using the Feynman rules in TFD J T n the 
thermal doublet notation the Feynman rule for the four-point scalar vertex is given 
by a vector like coupling constant which consists of the coupling constants for the 
non-tilde and the tilde fields. Thus the following factor is assigned to the cfr 4 vertex, 

A a = A C \ ) . (4-3) 



ki 

— ► 



/ \ 
-> i — ► \ — ► 



P \ k 2 Jp 



— ► 

p-k r k 2 

Fig. 1. 2- loop thermal self-energy in \(f> 4 interacting model 
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In Fig. [T] we show the lowest order Feynman diagram which contributes to the 
time evolution equation for the neutral scalar field. According to the thermal prop- 
agator (|2-19|) and the coupling constant (|4-3p . the self-energy is given by 

iUjj-" 12 (t Zl — t Z2 ;p) 

~ 3! J (27r) 2 (27r) 2 1 ' ' ^ 1 L Z2> K '2)- L, \ L zi h Z2iH)i 

(4-4) 

with q = p — k\ — ki- Substituting the explicit expression for the thermal propagator 
(HHH, we obtain 

ilTg 172 \t zi —t Z2 ;p) 

f rf2fcl d2k2 1 C4^ 

fk\ ,h,a {pZ2 ) fk2,i2,a(tz2 ) /?,i3;1 (^2 ) — /fcl ,M ,6 (^2 ) /fc2 ,12 ,6 (^2 ) /q,«3 ,6 (^2 ) 
fki,ii,a(tz 2 ) fk2,t2,a{tz2 ) fq,iz,a {t Z2 ) — /fei ,ii,b(^Z2 ) /fa ,i2,6 (^2 ) fq,i3,b V-%2 ) 

+ 0(t Z2 — t Zi ) e - i (- E fcl,n+' B fc2.*2+' E '3,»3)(^l _ * z 2) 

fki,ii,a(t zi )fk 2 ,i 2 ,a(tzi)fq,i3,a{tzi ) ~~ fki,ii,a(tzi )fk 2 ,i 2 ,a ifzj )fq,is,a (tzi ) 
fki,i\,bV'Zx) fk2,i2,bV'Z\) fq,i3,bV'Z\ ) — fki,ix,b\pz\ )/fc2,i2,b(^i )fq,i3,bV'Z\ ) 

where 

/ 5l l,a(t) = A,2,a(t) = 1 + n,(t), (4-6) 

/<Z,1,&(*) = 1 + Wg(t), /g,2,6(*) = Wg(t). 

We notice that unperturbed Bogoliubov parameters appear in the internal lines. 
The time dependent part in Eq. (|4-5|) has the following form, 

V(t - = 0(±(t - t'))e ±W( ^ < ' ) . (4-7) 
It is rewritten in a Fourier integral form, 

y (t _ t ') = _i / 1 e - ip0 (*-O. (4 . 8 ) 

Above po integral is simplified in the on-shell approximation. We set po on the 
denominator in Eq. fj4-8[) to the on-shell value, po = Up. 1 ' Thus we can perform the 
Fourier integration and obtain 

V(t - = -i5(t - -j— ^. (4-9) 

w„ =F VF =F 
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It is decomposed into the real and the imaginary parts. 

1 



V(t - t') = -i5(t -OP 



±2m5 (Lo p -W) 



(4-10) 



where P denotes the principal part of l/(w p =F W). Under the on-shell approximation 
the 2-loop thermal self-energy (|4-5p simplifies to 



i£™(t-t';p) 



n=i*2=i*3=i 



d 2 k ± d 2 k 2 



n) 2 (2tt) 2 8u kl uj k2 uj q 



1 



w p + -Efci.u + E k2l i 2 + E q ,i 3 + ie 

fki,ii,a(t )fk2,i2,a{t )fq,i3,ci(^ ) ~ fki,ii,b{t )fk2,i2,b{^ )fq,i3,b{t ) 

fki,ii,a\P ) fk2,i2,a(t' )fq,t3,a(^ ) ~ fk\,ii,b{t ) fk2,i2,b{t )fq,i3,b{t ) 
1 



Co',, 



E, 



<?,«3 



It 



fki,i 1 ,a(t)fk2,i2,a(t)fq,i3,a(t) fki , ii, a (^)/fc2,i2, a,(t)fq,i3, a {t) 
fk 1 ,i 1 ,b{t)fk2M,b{ t )fq,i3,b{ t ) -fk 1 ,i 1 ,b( t )fk2,i2,b( t )fq,i3,b( t ) 



(4-11) 



The perturbed propagator, Dh, is necessary to evaluate the right-hand side in 
Eqs. f|3-28j) and ([3v29j). Here we drop the higher order corrections and use the unper- 
turbed propagator, Dq, instead of the perturbed one, Dh- Substituting the thermal 
self-energy (|4-5p into the time evolution equations for the Bogoliubov parameter, 
p-28p and (|3-29p and replacing the perturbed propagator with the unperturbed one, 
we derive a equation with the structure of the Boltzmann equation for the A</> in- 
teraction model. We will call it Boltzmann equation, 



^2 2 2 2 



dt* 



d 2 ki d 2 k 2 



1 



il=li2=li3=l 

:cos{(-o;p + E kl , h + E k2ji2 + E q ^)(t x - t s )} 

\ + nH : p(tsY) fki,i\,a{ts) fk2,i2,a.{^s) fq,i3,a{t 
nH,p(ts) fki,ii,b(ts) fk 2 ,i2,b(ts) fq,i 3 ,b(ts) 



(2tt) 2 (27r) 2 16u} p oj kl uj k2 uj q 



(4-12) 



It should be noticed that both Eqs. (|3-28p and f|3-29[) provide the same expression 
(|4-12p . In the on-shell approximation the principal part in Eq. (|4-10p is canceled out 
from the time evolution equation. Thus the Boltzmann equation (|4T2p reduces to 



v p {t x ) 



3! 



EEE 

tl=l 12=1 13=1 



d 2 ki d 2 k 2 



IT 



(2tt) 2 (2tt) 2 8ujpU) kl u k2 uj q 
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x5(u p — Efa^ — Ek 2j i 2 — E qj i 3 ) 

y l + H/f i p(t a ;))/fc li j lja (t 3 ;)// C2i j 2ja (t X )/ 1 j ) j 3!a (i a ;) 

— nH,p{tx)fk 1 ,i 1 ,b(t x )fk2,i 2 ,b(tx)fq,i 3 ,b{tx) • (4"13) 

The right-hand side of this equation contains the delta function for the energy con- 
servation and the statistical factors. It corresponds to the collision term in the 
Boltzmann equation. 

§5. Numerical analysis of the Boltzmann equation 

The Boltzmann equation (|4T2p describes the time evolution of the particle num- 
ber distribution. We numerically solve it starting from the Bose distribution with 
the temperature, T, and the mass, mo, 

n ' it = 0) = ^Wrk^[ - (54> 

It is assumed that the scalar mass suddenly changes from mo to m. Then the 
state is no longer in equilibrium. Below we set parameters to mo = 5 x 10~ 2 /i, 
m = 4 x 10 -2 ^, T = n and A = fx, with an arbitrary mass scale, \i. It is expected 
that the distribution function, n p (t), approaches to the Bose distribution, n/ iP , with 
temperature, Tf. The final state temperature can be estimated at Tf = 0.988/j, from 
the energy conservation law. 

First we numerically solve the Boltzmann equation (|4- 13|) under the on-shell 
approximation. It is the differential equation in terms of the time variable. We 
employ the fourth order Runge-Kutta algorithm to solve it. In each steps of Runge- 
Kutta algorithm the momentum integral is performed by the second order Simpson 
integral. The radial part of the momentum integral is cut off at the scale, A = 20/i. 
In order to reduce the numerical error the integral interval is divided into two sub 
intervals, [0, 0.3yl] and [0.371, A]. The Simpson's rule is applied to each sub interval. 
The scattering with the particle in the final equilibrium state is introduced through 
the self-energy. 

For the numerical analysis we define the deviation from the final equilibrium 
distribution with the temperature, Tf, by 

5n p (t) = n p (t) — nf >p . (5'2) 

We expand the distribution function f(t) in Eq. (|4-13p in terms of Sn, 

W)W = 4 a(b) + 0(5n(t)), (5-3) 

where i = 1, 2 and 

fLa = fq,2,b = n f,1> ( 54 ) 



16 



Y. Mizutani and T. Inagaki 



In our setup the initial state temperature, T = jx, is close to the expected final 
state temperature, Tj = 0.988/x. Then we keep only the leading order 5n expansion 
in numerical calculations. The Bogoliubov parameters is fixed to nj^ for the internal 
lines in the self-energy. We solve the time evolution of the one for the external line, 
n H . 



0.00 ^□□□□□□□□□ggH0a0aOS00OOHH02SeSHiaE]00H2fflE]E]HH- 



-0.02 



o- -0.04 



-0.06 
-0.08 
-0.10 







3 



Fig. 2. Behavior of 8n v {t) in the on-shell approximation (|4-13[) . The triangle and square points 
show 8n p (t — 0) and 5n p (t = 1.0 x 10 2 /i -1 ), respectively. 




Fig. 3. Behavior of n p (t) in the on-shell approximation (|4-13p . The triangle and square points 
show n p (t — 0) and n p (t = 1.0 x 10 2 /i _1 ), respectively. The solid line represents the Bose 
distribution at T = Tf. 



Evaluating the Boltzmann equation (|4T3p in above assumptions, we obtain the 
time evolution of the distribution function. In Fig. [2] the behavior of 5n p (t) is shown 
as a function of the momentum p. At t = the particle number density with a lower 
momentum is much smaller than the estimated final one. We also draw the behavior 
of the particle number distribution, n p (t), in Fig. [3l As is expected, the particle 
number distribution approaches the equilibrium distribution with the temperature, 

Tf- 
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To evaluate the off-shell mode contribution we rewrite Eq. (|4T2p as 



dt s v{t s ) 



2 2 2 2 

3]" XT XT 



d 2 k\ d 2 k2 



il=li2=l «3=1 

sin{(— ujp + E kli u + -Efc 2 ,i2 + E g: is){t x — t s )} 
—u p + E kl ^i + E k2 ^ 2 + E q ^ 

(1 + nH,p(t s ))fk 1 ,il,a{ts)fk 2 ,i2,a(t s )fq,i3,a(ts) 
~ nH, P {ts)fk 1 ,il,b{ts)fk 2 ,i2,b{ts)fq,i3 ) b{ts 



(27T) 2 (27r) 2 l6u p U kl U k2 U q 



(5-6) 



The Boltzmann equation (|4T2p is reproduced by differentiating both the sides of 
this equation in terms of t x . This equation is satisfied for an arbitrary t x . Thus we 
obtain the differential equation, 



, 2 2 2 2 



d 2 k x d 2 k 2 



1 



(2tt) 2 (2tt) 2 16u p u kl u k2 u q 



11=1 42=1 13=1 ' 

sin{(— Up + E kl)il + -Efc 2 ,j 2 + E q: i 3 )(t x — t s )} 
+ E kl ^ 1 + E k2 i 2 + Eq^ s 

1 ,i 1 ,a{ts) fk 2 ,i 2 ,a(t s ) fg : i 3 >a (i s ) 
■ nH,p(t s ) fki,il,b{ts) fk 2 ,i2,b(t s ) fq,i3,b{ts) ■ 



(5-7) 



The time integral is dropped. Therefore the same numerical algorithm for the on- 
shell approximation can be applied to Eq. ()5-7j) . The time variable, t x , in Eq. (|5-7p 
represents the time interval from the initial moment. Here we set t x = 1.0 x 10 3 // -1 . 
At the limit, t x — > oo, the off-shell mode contribution disappears. 




-0.04 



0.04 



6cu/ju 

Fig. 4. Behavior of St{5oS) for t x -t s = 1000^~\ 

The integral kernel in Eq. (|5-7p is proportional to 

sh\{5u{t x -t s )} 



S t (5u) 



5u 



(5- 
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As is shown in Fig. UJ it has a peak at the on-shell limit, 5u — > 0, and frequently 
oscillates for a larger 6u). Thus it is enough to evaluate the integral near the on-shell 
limit. We restrict the integral interval, —12tt < 5u(t x — t s ) < 12ir. 



0.00 □□□□□□□□□□□gBHfiQaBBSOOBHHHHHHQQHZSZSHHHHHH- 
A K 

A 

-0.02 



a. -0.04 



-0.06 

-0.08 
-0.10 



3 

p/jl/ 



Fig. 5. Behavior of Sn p (t) with the off-shell contribution (|5-7|) . The triangle and square points 
show Sn p (t — 0) and Sn p (t = f .0 x 10 2 /i _1 ), respectively. 
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1 2 3 4 5 6 

Fig. 6. Behavior of n v (t) with the off-shell contribution ()5-7|) . The triangle and square points show 



n p (t = 0) and n p (t 
at T = Tf. 



1.0 x 10 /i ), respectively. The solid line represents the Bose distribution 



We numerically evaluate the Boltzmann equation (|5-7|) in above approximation 
and obtain the time evolution of the distribution function with the off-shell mode 
contribution. The behavior of 5n p (t) and n p (t) is illustrated as a function of the 
momentum p in Figs. [S] and El respectively. Hence, a similar behavior is observed for 
Eqs. ()4-13p and (j5-7|) . It seems to be difficult to distinguish the contribution from 
the off-shell mode in these figures. 

To evaluate the contribution from the off-shell mode we calculate the relaxation 
time, T p . According to the linear response theory, we define it by 

5n n 



5ri, 



(5-9) 
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It can be calculated by solving the Boltzmann equations (|4-13p and ([5 
600 




Fig. 7. The relaxation time, r p , in the on-shell approximation. The triangle and square points 
show r(t = 0) and r(t = 1.0 x lO 2 /^ -1 ), respectively. 
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r 1 300 
200 
100 


1 2 3 4 5 6 

Fig. 8. The relaxation time, t p , with the off-shell contribution. The triangle and square points 
show r(t = 0) and r(t = 1.0 x 10 2 /i _1 ), respectively. 

We plot the behavior of the relaxation time, Tp , ELS cl function of the momentum, 
p, in Fig. [7] under the on-shell approximation. The off-shell contribution is included 
in Fig. [8l The relaxation time is almost constant with respect to the time variable, 
t. A longer relaxation time is observed for a higher momentum mode. As is shown 
in Figs. [3] and [6] the distribution function n(t) is close to the final equilibrium state 
at t = 1.0 x lO 2 ^ -1 . It is consistent with the behavior of the relaxation time. For 
p < 2.0/U the relaxation time is smaller than 1.0 x 10 2 ^ -1 . The particle number 
distribution for a higher momentum mode is close to the final equilibrium state at 
t = and slowly approaches the final state. 

A small but stable discrepancy is observed for some points around p ~ in 
Figs. [7]and [3 It corresponds to the numerical error which is produced due to a finite 
step size in the Simpson integral. 
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Fig. 9. Behavior of the relaxation time, t p . The triangle and square points show the relaxation 
time without and with the off-shell contribution, respectively. 



In Fig. [9] both of the results are drawn in the same figure. The off-shell mode 
contribution does not modify the relaxation time for a lower momentum. For p > 3fi 
we observe a shorter relaxation time in the on-shell approximation. Since the off- 
shell contribution destabilizes the particle distribution near the equilibrium state, 
it slightly increases the relaxation time. A different contribution from the off-shell 
mode is expected for a system far from the equilibrium state. Because of a non- 
negligible numerical error it is difficult to apply our analysis to such a system. Some 
improvements of the numerical algorithm are necessary to extend our analysis for a 
general case. 

§6. Conclusion 

We have investigated the time evolution of the distribution function for a rela- 
tivistic neutral scalar field with a Ac/) 4 interaction. The NETFD is applied to the SD 
equation for the scalar field propagator. Calculating the 2-loop thermal self-energy 
and inserting it into the SD equation, we have derived the time evolution equation 
for the thermal Bogoliubov parameter. The equation has the same structure with 
the Boltzmann equation. Therefore the Boltzmann equation is obtained from the 
SD equation at the equal time limit. 

The Boltzmann equation consists of the collision terms which depend on the 
thermal Bogoliubov parameter and an oscillating coefficient with respect to the time 
variable. The oscillating coefficient has a peak for the on-shell case. It reduces to the 
delta function which shows the energy conservation between the collision particles 
in the on-shell approximation. Solving the obtained Boltzmann equation, we have 
evaluated the time evolution of the thermal Bogoliubov parameter. We suppose 
that the scalar field loses 20% of mass suddenly at t = 0. It is observed that the 
thermal Bogoliubov parameter approaches an equilibrium state. The relaxation time 
monotonically increases as a function of the momentum. In our setup the off-shell 
mode has only a small contribution to the particle number density. The off-shell mode 
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tends to suppress the relaxation near the equilibrium state and slightly increases the 
relaxation time for a higher momentum region. 

We have derived the Boltzmann equation from the SD equation. For a non- 
relativistic scalar field the Boltzmann equation has been derived as the self-consistent 
renormalization condition based on the canonical quantization by H. Umezawa and 
Y. Yamanaka.E'M In this case it has been shown that the thermal Bogoliubov 
parameter which satisfies the self-consistent renormalization condition coincides with 
the observed particle number densityP'^'ED We would like to show the correspon- 
dence between the thermal Bogoliubov parameter and the particle number density 
in our approach. For this purpose it is necessary to construct our procedure based 
on the canonical quantization. 

In the derivation of the time evolution equation the second derivative of the 
Bogoliubov parameter is canceled out at the equal time limit. Thus we obtain 
the Markovian equation. In the calculation of the self-energy part we impose the 
internal scalar propagator to be in the final equilibrium state. We should improve 
these assumptions to generalize the procedure. It is also interesting to apply our 
analysis to the relativistic Dirac field. These will be the subject of a forthcoming 
paper. 



The authors would like to thank Y. Yamanaka, H. Matsumoto and Y. Naka- 
mura for valuable discussions. Discussions during the YITP workshop on "Thermal 
Quantum Field Theories and Their Applications 2010" were useful to complete this 
work. 



The relationship between the perturbed and the unperturbed propagator is given 
by the SD equations (|3-ip and (|3-2p . Here we evaluate the diagonal elements of these 
equations. Applying the Bogoliubov matrix, i? _1 (n#) and B(nn) to Eq. (|3-ip on the 
left and the right sides, respectively, we obtain conditions for each diagonal element 
of the perturbed propagator. 
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Appendix A 

diagonal elements of the SD equation 





+ J dt Zl dt Z2 (D$ R (t x -t Zl ) + Dl 2 A (t x - t Zl )) 

X iS R (t zl - t Z2 )(D l l R (t Z2 - ty) + D% A {t Z2 - ty)) , 



(A-l) 



^H,R(^X — ty) + D HA (t X — ty) — D QR (t X — ty) + D HA (t X 




+ J dt zl dt Z2 (Dl 2 R {t x -t zi ) + Dl] A {t x - t Zl )) 

X iZ A (t Zl - t Z2 )(D 2 l R (t Z2 - ty) + D X l A (t Z2 - ty)). 
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We suppose that the retarded and the advanced parts of these equations separate 
and impose conditions, 

D% R (t X - ty) = Dl] R (t x - ty) (A3) 

+ j dt Zl dt Z2 (^D^ R (t x - t zi ) + Dl 2 A (t x - t Zl ))iU R (t zl - t Z2 )D^ R (t Z2 - t y ), 
D 2 i A (t* ~ t y ) = D 22 A (t x - t y ) (A-4) 
+ J dt Zl dt Z2 (D$ R (t x - t Zl ) + Dl%t x - t Zl j)iE R {t Zl - t Z2 )D% A (t Z2 - t y ), 

D% R (t X ~ ty) = Dl 2 R {t x - ty) (A-5) 

+ J dt Zl dt Z2 (Dl 2 R (t x - t Zl ) + Dl] A {t x - t Zl fjiE A {t Zl - t Z2 )D 2 l R (t Z2 - t y ), 

D% A (t X - ty) = Dl] A {t X - ty) (A-6) 

+ J dt Zl dt Z2 (D 22 R {t x - t Zl ) + Dl] A {t x - t Zl fjiE A {t Zl - t Z2 )D]} A (t Z2 - t y ). 

From the SD equation (|3-2p we find other conditions for the perturbed propagator. 
DH,n(tx — t y ) + D 2 x A {t x — t y ) = D 0R (t x — t y ) + D^ A (t x — ty) 
+ J dt Zl dt Z2 (D% R (t x -t Zl ) + D 2 l A (t x - t Zx )) 
x iU R (t Zl - t Z2 )(Dl] R (t Z2 - t y ) + D 22 A (t Z2 - t y )) , (A-7) 

^H,R^ — ty) + D HA {t X — ty) = DQ 2 R (t X — ty) + D HA (t X — ty) 

+ J dt zl dt Z2 (D 2 £ R {t x -t Zl ) + D% iA (t x - ttS) 

x iZ A (t Zl - t Z2 )[D 22 R {t Z2 - t y ) + Dl] A {t Z2 - t y )) . (A-8) 

In a similar manner with Eqs. (|A-3P - (|A-6p we divide Eqs. (|A-7j) and ()A-8p . 

D 1 i R {t x -t y ) = Dl] R {t x -t y ) (A9) 
+ J dt Zl dt Z2 D]} R (t x - t Zl )iU R (t Zl - t Z2 )(^D^ R (t Z2 - ty) + Dl 2 A (t Z2 - t y )j, 

D% A {t x ~ t y ) = D 22 A (t x - t y ) (A-10) 
+ J dt Zl dt Z2 D 2 x A (t x - t Zl )iE R (t Zl - t Z2 )(^D^ R (t Z2 - t y ) + D^ A (t Z2 - t y j) , 

D 2 H,R(t x - t y ) = D 22 R (t x - ty) (A-ll) 
+ J dt Zl dt Z2 D 2 x R (t x — t Zl )iU A (t Zl — t Z2 )^p^ R {t Z2 — t y ) + D^ A (t Z2 — ty)) , 

D l i A {t X - ty) = Dl) A {t X - ty) (A-12) 

+ J dt Zl dt Z2 D]j A (t x — t Zl )iU A (t Zl — t Z2 )(^DQ 2 R (t Z2 — t y ) + D^ A (t Z2 — t y ) 
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The self-energy, S r< a has retarded and advanced time dependence, En(t — t') = 
for t < t' and 27a — t') = for t > t'. The time dependence for D R and Da is given 
in Eqs. (|2-2U|) - (|2-23p , Thus the quantum corrections for the perturbed propagators 
(|A-3p - (|A-6p and ()A-9[) - ()A- 12|) disappear at the equal time limit. As a result the 
perturbed propagators at the equal time limit are derived from the unperturbed 
propagators (l2^ - (l2 : 23l) . 

lim D% R (t x - t y ) = lim D% A (t x - t y ) = (A-13) 

lim D% R {t x - t y ) = lim D% A (t x - t y ) = (A-14) 

t X rty t X rty QuJp 

Substituting the expressions (|2-20p - (l2-23p for the unperturbed propagators, we 
obtain the first derivative of the perturbed propagators at the equal time limit, 

lim B tx D™ R {t x - t y ) = lim d ty D% A (t x - t y ) = % - + lim 5(t x - t y ), (A-15) 

t x rty t x rty Zi t x rty 

lim d tx D% A (t x - ty) = lim d ty D x ^ R (t x -t y )= % -- lim 5(t x - t y ). (A46) 

tx TLy tx tty £ tx rty 
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